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Abstract. We study the time-dependent Gross-Pitaevskii equation describ- 
ing Bose— Einstein condensation of trapped dipolar quantum gases. Existence 
and uniqueness as well as the possible blow-up of solutions are studied. More- 
over, we discuss the problem of dimension-reduction for this nonlinear and 
nonlocal Schrodinger equation. 



1. Introduction 

The success of atomic Bose-Einsteiii condensation has stimulated great interest 
in the properties of trapped quantum gases. Recent developments in the manip- 
ulation of such ultra-cold atoms have paved the way towards Bose-Einstein con- 
densation in atomic gases where dipole-dipole interactions between the particles 
are important. In [21] Yi and You were the first to introduce a pseudo-potential 
appropriate to describe such systems in which particles interact via short-range 
repulsive forces and long-range (partly attractive) dipolar forces. Describing the 
corresponding Bose-Einstein condensates within the realm of the Gross-Pitaevskii 
(mean field) approximation, one is led to the following nonlinear Schrodinger equa- 
tion for the macroscopic wave function of the condensate [13, IG, 17, 22, 23]: 

(1.1) ihdti^ + — Ai/' = y(a;)V' + g|V'|V + f^^(i^*lV'nV', xeM^^>0, 

where \g\ = ATifi^N\a\/m. Here, we denote by N the total number of particles 
within the condensate, whereas m denotes the mass of an individual particle and a 
its corresponding scattering length (which can be experimentally tuned to be either 
positive or negative). The wave function ip is then normalized, such that [jV'llia = 1. 
The potential V{x), for x — {xi,X2,X7i) E M'^, describes the electromagnetic trap 
for the condensate and is usually chosen to be a harmonic confinement, i.e. 

(1.2) V{x) - I {cj'^xl + ^Ixl + u;lxl) . 

The real- valued constants wi , 0^2 , W3 then represent the corresponding trap fre- 
quency in each spatial direction. Finally, d denotes the dipole moment (in Gaussian 
units) and 

„x N 1 ~ 3cos^ 9 

(1.3) K{x) = — , 
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where 9 stands for the angle between a: G M"^ and the dipole axis n e M'^, with 
\n\ = 1. In other words 9 is defined via 

X ■ n 



cos 9 



X 



To our knowledge a rigorous mathematical study of (1.1) has not been given so far. 
Note that the interaction kernel (1.3) is indeed highly singular and it is therefore 
not clear a priori if the corresponding convolution operator is well defined. In 
particular, not even the existence and uniqueness of solutions to (1.1) has been 
established yet, and it will be one of the main tasks of this work to do so. Moreover, 
we shall be interested in the mathematical problem of dimension reduction. Quasi 
two-dimensional (pancake shaped) or even quasi one-dimensional (cigar shaped) 
Bose-Einstein condensates can be obtained experimentally by appropriately tuning 
the trap- frequencies u)i,uj2,(^3- The mathematical derivation of effective models in 
lower dimensions via the corresponding scaling limits is therefore of great practical 
importance. 

For the mathematical analysis it is more convenient to rescale (1.1) in dimen- 
sionless form (see e.g. [3]), to arrive at the following model 

(1.4) idtij + ^Ai' = V{x)i' + Xi\tlj\^tl; + X2{K*\ip\^)'4j, xeR^, 



where Ai = AiraN/aQ and A2 = (P / (hujoaQ) . Here, we denote by ao = yh/nwjQ the 
ground state length of a harmonic oscillator corresponding to loq = \ {oJi -I-W2 +W3). 
In the upcoming analysis Ai, A2 will simply be assumed to be two given, real- valued 
parameters. 

The paper is then organized as follows: We start by collecting several properties 
of the interaction kernel K in Section 2. After that, we set up a local in time 
existence theory in Section 3, before proving several different global existence results 
in Section 4. Section 5 is devoted to the question of finite time blow-up of solutions 
and also includes an additional global in time existence result for rather particular 
circumstances. Finally, we shall study the dimensional reduction of our model in 
Section 6. 

2. Some properties of the dipole kernel 

To simplify notations, we shall from now on assume, without restriction of gen- 
erality, that n = (0, 0, 1). The dipole-interaction kernel K then reads: 

(2.1) Kix) = 

Even though this kernel is highly singular (like l/jccp), it defines a rather smooth 
operator. The technical reason for this is that the average of K vanishes on spheres. 

Lemma 2.1. The operator K, : u ^-^ K can he extended as a continuous operator 
on LP(M3) for alll<p <oo. 

2 

Proof. We notice that O : x i-^ 1 — Scos^^ = 1 — 3j^ is Lipschitzean on S^^ 
homogeneous of degree zero, and 

for 2; G R'^. We also note that the average of H, on spheres vanishes: 



The lemma then follows from the Calderon-Zygmund Theorem (see e.g. [18]). □ 
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Remark 2.2. The fact that the average of K vanishes on §^ imphes that the 
dipole-nonhnearity vanishes when apphed to radially symmetric wave functions 
■0 = -^da;!). In this case (1.4) simplifies to the classical (cubic nonlinear) Gross- 
Pitaevskii equation. 

The continuity at the L^-level can also be seen by computing the Fourier trans- 
form of K, which turns out to be essentially bounded, i.e. K e L°°(IR'^). In the 
following we shall use the explicit formula for this Fourier transform several times, 
see also [13, 16]. 

Lemma 2.3. Define the Fourier transform on the Schwartz space as 
Tu{C) = u{0 = I e-"-^u{x)dx, u G 5(R^). 



Then the Fourier transform of K is given by 



|,3„.e-,.|(4-.).f(M_|-^^ 



where Q stands for the angle between ^ and the dipole axis n ~ (0, 0, 1). 
Proof. We use the decomposition of e~*^'^ into spherical harmonics 

oo 

^ ^ (_j)™(2m + l)j„ (|:e||^|) P™ (cos </.) , 

where jm and P,„ stand for the spherical Bessel function and Legendre polynomial of 
order m, respectively, and (j> denotes the angle between x and ^, i.e. a;-^ = \x\ |^| cos (j) 
(see e.g. [1]). Denoting by uj and uj' the corresponding spherical coordinates of x 
and ^, respectively, we shall further expand Pm via 



. m 

A'jr 



where Yim are the spherical harmonics. Note that, since 

P2{x) = i (3x2 _ ^ 

we have the following identity 



where we recall that 6 is the angle between x and the dipole axis n. 
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Putting all of this together we can compute the Fourier transform of K via 
K{0 V(-*r(2m+l)>„(|x||e|)P™(cos0)X(x)dx 

/ A ^ POO -| 

/ {cos 4,)Y2o{uj)dS{uj) 

2\/^EHr(2™+i) / j™wei)-x 

m 

^ y,:„ (a;) r,™(a;')>2o(^) dS{uj) 



m=0 

47r 



2m + 1 

/ 47r , ,9 , fir 47r , 

= - 2a/ — X 5 / ]2{r\£.\) - X —Y2a{u'), 
V o Jo r 

where in the last equality we have used the fact that the {Y^m} form an orthonormal 
basis of i^(§^). We thus have, for C 7^ 0, 

^ /ItT f°° dr f°° dr 

m) = SttJ— y2o(c^') / J2(r|^|)- = 87rP2 (cos 6) / j^M- 
V Jo T Jo ''' 

f°° dr 

= 47r(3cos2e- 1) / j2{r) — . 

Jo "f" 

It remains to compute the radial integral in this formula. In view of the identity 

_^ f h{r) \ ^ .72 (r) 
dr \ r J r ' 

we have 

J2[r) — = hm / j2[r) — = hm 



Jj^ r R^o R 

Recalling that 

sini? cosi? 

MR) JT' 

we thus find j2ir)/rdr — 1/3, which finishes the proof. □ 

3. Local in time existence 

In view of the forthcoming dimension reduction analysis in Section 6, we allow 
the spatial dimension to be smaller than three and consider, instead of (1-4), the 
following initial value problem 

^g^^ UdtiJ+^Ai,^Vd{x)^ + XM^ij + X2{Kd*\ij\^)ij, xGM^ 

where 1 ^ d ^ 3. We assume Ai, A2 £ M and the potential Vd to be quadratic in d 
dimensions, i.e. 

1 

(3.2) Vdix) = -Y,^]xl u;,eR. 

3 = 1 

The interaction kernel Kd will not be specified in detail (except that K^, = K). 
Rather, we shall only assume that the operator ICd '■ u 1— > Kd * m is bounded on 
L^(M'^), a property already known to be satisfied in d = 3, see Lemma 2.1. In 
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d = 1,2 we shall check that this holds true after we succeeded in deriving the 
precise expressions of Kd via dimension reduction. 

The following two important physical quantities are formally conserved by the 
time-evolution associated to (3.1): 

(3.3) Mass: M =||V^(t)||i2. 

(3.4) Enci-gy: E=hvm\\h+ [ Vd{x)m,x)\^dx + ^\m)\\l^ 

+ ^ [ {Kd*\^\'){t,x)m,x)\'dx. 

Note that for the conservation of energy (which can be derived formally by multi- 
plying (3.1) by dttp, taking real parts and integrating in x), we use the fact that K 
is even. These two quantities naturally lead to the introduction of an energy space 
associated with the linear case Ai = A2 = 0: 

S = {« e L^iR") ; \\u\\l := + \\\/u\\l, + \\xu\\l, < ^} . 

Let us ignore the dipolc nonlincarity for the moment: A2 = 0. S is a natural space 
to study (3.1), where several results concerning the Cauchy problem arc available 
(see e.g. [!)]). Moreover, it is known that working only in L^(R^) is not sufhcient to 
prove (local) well-posedness; see e.g. [7] and references therein. Also note that the 
Cauchy problem is ill-posed in H^{1S.'^), due to the rotation of phase space induced 
by the harmonic oscillator (see [S] for both linear and nonlinear cases). However, 
including the dipole-nonlinearity in this setting is not completely straightforward as 
we can not allow for any (weak) derivative to be put on K^, since the corresponding 
convolution operator is no longer well defined. The only possible way to circumvent 
this problem seems to be the use of Strichartz estimates. To this end, we shall 
introduce in the following subsection the main technical tools needed later on in 
the local existence proof. 

3.1. Technical preliminaries. Let 1 ^ d ^ 3, and Vd be quadratic as in (3.2). 
Introduce the group 

U{t) = e""^, where H = -^A + Vd, 

which generates the time-evolution for the linear problem. We first remark that 
in view of Mchlcr's formula (see e.g. [10]), the group [/(•) is not only bounded on 
L^{M.'^), but also enjoys dispersive properties for small time. More precisely it holds 

C 

(3.5) |lC/(i)^|lz,.^(Rd) s$ |^P7jI!^IIli(r^), for \t\ ^ 5, 

for some (5 > 0. Note that for an harmonic potential as in (3.2), S is necessarily 
finite, since H has eigenvalues. 

Remark 3.1. This dispersive estimate (3.5) turns out to be valid for external 
potentials which are not exactly quadratic (nor even positive): it suffices to have 
Vd e C°°(R'^;R) with d^'Vd e L°°(M.'^) as soon as \a\ ^ 2 for the above estimate to 
remain valid. In that case, the proofs we present below remain valid, see e.g. [6] 
and references therein. 

An important consequence of the dispersive estimate (3.5), and of the fact that 
[/(•) is unitary on i^, is the existence of Strichartz estimates, which require the 
following definition of admissible index pairs. 
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Definition 3.2. A pair {q, r) is admissible if 2 ^ r ^ (resp. 2 ^ r ^ oo if 
(i=l, 2^r<(X)if(i = 2) and 

q \2 r 

Following [11, 20, 15], we have: 

Lemma 3.3 (Strichartz estimates). Let {q,r), (qi,ri) and (92, ^'2) be admissible 
pairs. Let I be some finite time interval containing the origin, € /. 

1. There exists Cr = C{r,I), such that for any (p £ L^(]R'^) it holds 

(3.6) \\U{-MLHI;Ln'^'^rML^- 

2. There exists Cri,r2 = C!{ri,r2, 1), such that for any F G L'^2(/; ^^''2) a holds 



(3.7) 



/ U{t - s)F{s)ds 

Jin{s<t} 



L9i(/;L'-i) 



As underscored above, taking for (p an eigenfunction of the (anisotropic) har- 
monic oscillator shows that in general the constants Cr and Cri,r2 do depend on 
the length of the time interval /, unless all the w/s in (3.2) are zero. 

In the following we shall use two special admissible index pairs to apply the 
above lemma, namely 



(3.8) (g,r) = (oo,2) and (g,r)=(^-,4 

The idea for the proofs of local in time existence and uniqueness is then to apply 
a fixed point argument on the Duhamel's formula associated to (3.1). This step 
turns out to follow exactly the same lines as in the proof of local existence without 
dipole, A2 = 0. To this end Strichartz estimates based on (3.8) will be invoked 
several times. 

3.2. Constructing a local solution in the energy space. With the above 
technical tools in hand we can now state the first main result of this work. 

Proposition 3.4. Let 1 ^ d ^3, Vd be quadratic, Ai, A2 G K, and ip gT^. Assume 
that the operator JCd ■ u 1— *■ Kd * u is bounded on L^(R'^). Then there exists T, 
depending only on (upper bounds for) \\^p\\y,, such that (3.1) has a unique solution 
ip G Xt, where 



[^gC{[Q,T]-^) ■ iJ,Vi;,xi^GC{[0,T]■L\W'))C^L^/\[Q,T]■L\W'))] 



Moreover, the mass M and the energy E, defined in (3.3), (3.4), are conserved for 

te [o,r]. 

The space L^/'^ ([0,r];L* (M'')) is here to ensure uniqueness. The appearance 
of this space will become clear during the course of the proof. The main technical 
remark is that since JCd is continuous on L'^{W^) by assumption, the nonlocal term 
in (3.1) can be estimated like the local cubic nonlinearity. 

Proof. The proof of the conservations of mass and energy is omitted here. We refer 
to [9] for a general argument, easy to adapt to the present case. Now, consider 
Duhamel's formula associated to (3.1) 

m^U(t)^-iXi [ Uit-s){\ij\^ij){s)ds~i\2 f U{t~s){{Kd*\^\^)iP){s)ds, 
Jo Jo 
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and denote ^{ip){t) = U{t)ip + Si{t) + S2{t), with 

Si{t)^-tXi f U{t-s){\ij\^^)is)ds, 







S2{t) = -1X2 I U{t - s) {{Kd * ^) {s)ds. 







Since at this stage, we are interested in a local result only, we may assume that 
T e]0, 1]. For 1^ e E, let i? = ||v3||s, and introduce 

Xt{R) ^{ijeXr ; || A?M|ioc,([o.T];L2) 4C2i?, 

and ||A^/'||LS/d([o,T];LJ) =^ 4C4i?, Vy4 e {Id, V,a::}}. 

The constants C2 and C4 arc those which appear in the first point of Lemma 3.3, 
with / = [0, 1]. Apart from them we shall in the following denote all the irrelevant 
constants by C, whose value may therefore change from one line to another. 

To prove Proposition 3.4, we first show that for T €]0,1] sufficiently small, $ 
leaves Xt{R) invariant. Then, up to demanding T to be even smaller, we show 
that $ is a contraction on Xt{R), for the topology of L*/''([0, T]; L^). As remarked 
in [14], Xt{R), equipped with this topology, is complete, and thus Proposition 3.4 
follows, since T will only depend on d and R. 

Step 1 (stability): Let V e XriR) and denote Lf.L'' = L«([0, T]; i''(M'')). The 
Strichartz estimates yield: 

II^WIIl^l^ \\U{-Ml^l^- + ||5i|U=oi. + \\S2\\l^l^ 

< C2R + 6*2,4 II \lp\'^i'\\^8/(8-d)^i^3 + 6*2,4 II {Kd * IV'P) V'lliS/cs-dj^i/a . 

Similarly, it holds 

II^WIIiB/d^, 5:C4i? + C4,4|||^PV'|L + 6*4,4 II (A'dHV'l')V'IL 

We shall only work on the latter estimate as the L^L^ estimate can handled exactly 
in the same way. Denote (g, r) = (S/d, 4) and remark the following identities 

3 13 8-(ill2.,, 8 
3.9 7 = - = - ; ^ = - = - + 1:' ^ithfc-— — . 

4 r r 8 q' q k A — a 

Holder's inequality thus yields 

l|l^l^^ll4/<«-'')LV3 < Ikili^.L^ll^lli'L'-- 

By assumption, the operator ICd ■ u Kd * |up is bounded on L^(R'^) and so we 
have (recall that r = 4): 

We note that this is exactly the same estimate as for the local cubic term, up to 
the norm of ICd- From Sobolev embedding (in space), we get 

II^IIl^l- ^ CII^IL.^, < CT'^'^Ml^h^ ^ CT'/^C2R. 

Therefore, 

mmL^^Lr^c,R+cT^/''R\ 

and if T e]0, 1] is sufficiently smaU, ||$(V')||l^L'- 46*4i?. 

Next, to estimate V<f>(7/;), we denote by [A, B] = AB — BA the usual commutator 
and compute 

[V, H] = [V, V] - Vy, [x, H] = A] = -V. 
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By assumption, |VF(x)| ^ C\x\, which shows that we wih indeed get a closed 
family of estimates for and For instance, 

V$(V')(i) = U{t)Vif - i\i [ U{t - s)V (IV-lV) {s)ds 



(3.10) -1X2 U{t-s)V{{Kd*\i^\^)i^)is)d. 



U{t-s) ($(V'(s))VF) ds. 
From the Strichartz estimates, we thus have 

I|c/(-)v^||l-l2 < C2R, 

and for the local cubic nonlinearity we get 



U{t- s)\7 {s)ds 



where we have used the same computations as above, on For the dipole- 

nonlinearity, we similarly obtain 

t 

U{t - s)V {{Kd * ^) {s)ds ^^^^ ^ C24 II V {{Kd * IVP) ^) IL^'i,. 

< C II {Kd * IVP) VVlL,'^.' + C II (Kd * VIVP) V'lL,'^.' 

< C\\Kd * |^n!i./2^./.||VV||L^L. + C\\Kd * vivpll 

(l/<J+l/fc)-l ;-r/2 

IIV-ILu- 

This is the same type of estimate as for the local cubic term. Finally, the last term 
in (3.10) is estimated by 

U{t~s){m,{s))\/V)ds ^^^^ 5:C2,2||$(V')VFL^i2 ^CT||a-$(V)|U^i2. 

All in all, we find: 

I|V$(V)||l~l2 C2R + CT^'^'R' + CT\\x(^m\L^L^.. 
Analogously we obtain (with slightly shorter computations) 

WxH^L^L^- < C2R + CT^'^R^ + CT\\W<^m\L^L^-. 

Summing up these two inequalities and taking T g]0, 1] sufficiently small, we arrive 
at 

I|V$(V)|1l-L^ + WxHniL^L- ^ 4C2i?, 

and using the Strichartz estimates again, wc also have 

I|V$(V)IIl^l^ C4i? + Cr2A-i?3 + CT||a;$(V)|U^.i., 
W^'^Wh^^L'- < CiR + CT'^'R' + CT||V$(V)|U5?L- 

Up to taking T even smaller, we therefore see that $ leaves Xt{R) stable. 
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Step 2 (contraction): For ipi and 'ip2 in Xt{R), we have, from Strichartz estimates 
and the above computations, 

||$(^i)-$(V'2)||l^L^ 

^ C4,4lll^iPV'i - \^2\^H\l^;l^' + CiA {Kd * IViP) i'l - {Kd * ^2\\^,-^., 
^ C II (iViP + 1^1 - MWl^;^^' + C II {Ka * (l^ip - |^2p)) MLiL^' 

+ c||(ifdHV'2n(V'i-V'2)||^^'^.- 

^ CT2Ai?2|lV'i - V'2||l^l^ + C lllViP - |V2nLa/,+iA-,-^. IIV^iIIl^l^ 

+ CT2/^i?2||V.i-V^2||L^L. 

Therefore, $ is a contraction provided that T is sufficiently small, and hence the 
resuh follows. □ 

4. Global in time existence results 

Having established the existence and uniqueness of solutions to (3.1) locally in 
time, we now turn our attention to global in time results. To this end we first note 
that the only obstruction for global existence in S is the possible unboundcdncss of 
VV' in L?{W'-). Indeed, in W'- with 1 ^ d ^ 3, the Gagliardo-Nircnberg inequality 
(see e.g. [!)]) yields 

\\u\\\.^C\\ut-/\\Vu\\i... 

Thus, the boundedness of JCd on L'^{W^) together with the Cauchy-Schwarz in- 
equality implies 

^^^^ ,_ {Kd*\u\^){x)\u{x)\^dx 



^C\\ut-/\\Vu\\i,, 



where we have used Gagliardo-Nirenberg inequalities. So we see that in view of the 
conservation of the mass, the (conserved) energy is the sum of four terms, three of 
which are bounded provided that Vij) remains bounded in L'^{W'-). Therefore, the 
fourth term is bounded, that is, xtp{t,x) is bounded in L^(R''). This shows that 
unless ||VV'(i)lli2 becomes unbounded, ||'0(i)||s is a continuous function of i ^ 0. 
This directly yields the following corollary: 

Corollary 4.1. Under the same assumptions as in Proposition 3.4, there exists a 
G U {00}, such that (3.1) has a unique maximal solution in 

[iP e C([0,T4;S]) ; ^, VV',.tV e C([0, [; ^^(R'^)) n Lf/;^([0, T,[; ^^(m'^))} , 

such that "014=0 = H is maximal in the sense that if < 00, then 

l|VV(i)llL2^^^+oo. 

Having in mind this result, we are now able to study the global in time existence 
of solutions to (3.1) depending on the spatial dimension. In the following we shall 
first treat the case d = 1, where the picture is much more concise, before moving on 
to the case d = 3 (the case d = 2, is similar to the one of three spatial dimensions 
and thus omitted for simplicity). Note that, in view of the dimensional reduction 
discussed in Section 6, the one-dimensional case is not purely academic. 
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4.1. Global existence for d = 1. When d = 1, the solution constructed in the 
previous section indeed exists for all time. More precisely we have: 

Corollary 4.2. Suppose that in Proposition 3.4, d = 1. Then the solution to (3.1) 
is global in time, i.e. = oo in Corollary 4-1- 

Proof. This result stems from Gagliardo-Nirenbcrg inequality, like its classical 
counterpart when A2 = 0. In view of Corollary 4.1, we just need an a priori 
estimate for |j VV'(^)||l2. Indeed the conservation of energy yields 



\\vm\\h^2E+\x,mt)\\h + \>^2 
< 2E+c\m)\\'U\^m\\L- 



{Ki*\i^\^) {t,x)m,x)\''dx 



where we have again used the estimates (4.1). Since M = ||V'(0lli2 is constant in 
time, this inequality directly shows that || V^/'(t)||i2 remains bounded and hence 
the result is proven. □ 



When d = I, a more general setting is possible, which allows us to work with 
initial data that are only in but not necessarily in E (corresponding to solutions 
with possibly infinite energy). The following result is an adaptation of the main 
result in [19]. 

Theorem 4.3. Let d ~ 1, Vi be quadratic, and Ai,A2 G K. Assume that the 
operator K-i : u t-^ Ki * u is bounded on L^(M). Then, for any ip G L^(]R), (3.1) has 
a unique solution 

^eC{R+;L^)nLl, {R+;L^ (M)) . 
Moreover its total mass M is independent oft ^ 0. 

Proof. Wc shall first prove that there exists a T, depending only on ||(y3||^2(jj), such 
that (3.1) has a unique solution 

VeC([0,T];L2)nL«([0,r];i4 (R)), 

To this end we resume the same scheme as in the proof of Proposition 3.4 above, 
but now R = ||iy9||L2 and 

rT(i?) = {7^eL°°([0,T];i2)nL8([0,r];i4) ; ||V^|Uc.([o,T];L2) ^ 2C2i?, 
and ||V'IIls([o,t];L-») 2C4R}. 

We simply notice that for d = 1, the identities (3.9) yield k = 8/3, and so Holder's 
inequality in time implies 

Following the same lines as in the previous paragraph, we sec that $ leaves Yt{R) 
stable, provided that T is sufficiently small. The proof of contraction is the same as 
in the previous paragraph, up to the modification of the estimate for HV'IIl^l'^- This 
yields a local in time existence and uniqueness result and since the existence time 
depends only on ||(y9||^2(ig), the conservation of the total mass M directly implies 
the global in time existence. □ 

In summary the Gross-Pitaevskii equation (3.1) is globally well-posed for = 1. 
The three dimensional case is more involved, though, as we shall see. 
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4.2. Global existence for d = 3 (stable regime). We can now turn to the 

physically most important case d = 3. 

Theorem 4.4. Under the same assumptions as in Proposition 3.4, suppose that in 
addition Ai ^ ^0. Then the solution is global in time, i.e. T* = oo. 

In the following, the situation where Ai ^ ^0 will be called the stable 

regime (note that A2 > corresponds to the actual physical situation). 

Proof. Wc first note that from Planchcrcl's formula for p = \^\'^ we get 

\\m\\U = \\pmh = j^\m\\l^- 

Then we simply use the conservation of the energy E, as defined in (3.4), to estimate 

wvmwh^^E- f V3{x)m,x)\^dx ~ x,\m)\\U 

-A2 / {Ki*\i'\-'){t,x)m,x)\''dx 

invoking the above given identity. Recalling the explicit formula for K computed 
in Lemma 2.3, we obtain 

By assumption, this quantity is non-negative and hence the a priori estimate 

holds true. The result then follows directly from Corollary 4.1. □ 

Having established this result, it is natural to ask what happens in the unstable 
regime Ai < ^A2. As will shall see, in general we cannot expect a global in time 
result there, since finite time blow-up of solutions (in the sense of Corollary 4.1) 
may occur. In particular, it is clear that if blow-up occurs, then the problem of 
dimension reduction ceases to make sense. The next section is devoted to the study 
of this problem. 

5. The unstable regime 

When d = 3 (or 2), the solution constructed in Proposition 3.4 need not remain 
in S for all time. This will be seen from using a general virial computation (see [9]) 
and following the approach by Zakharov [24] and Glassey [12]. 

5.1. Finite time blow-up. As a preliminary, we check that the energy E may be 
negative. 

Lemma 5.1. Let d — 3. Assume that 

Ai < yA2. 

There exists ip d Y, such that E < 0, where the energy E is as in Proposition 3.4-. 
Proof. We first rewrite the last term in the energy, thanks to Planchercl formula: 

/ {K^\^\^){^xmxrdx = -^ f m)\m?di. 
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where we have denoted p = l^pp. In view of Lemma 2.3, we infer 



{K*\^\'){x)\ip{x)\'dx 



1 



R3 



■ii 



1 \m)?dc 



The idea of the proof then consists in choosing so that p has httle mass on 
{^3 = 0} = • n = 0}. Indeed, using Plancherel formula again, we can rewrite the 
energy as 



(5.1) 



2(27r) 



2) and 5 e S{R) 



Introduce a parameter e > 0, and fix some functions / G 
independent of e. Set 

if{x) ^ e"-/'^f{xi,x2)g{ex3), 
for some constant a to be fixed later. We have 



p{x) = Mx)\' ^s"\fix^,X2)\'\gisx3r , -£"-^^^(6,6)6 



^ e 

where F and G denote the Fourier transforms of |/p and \g\'^ ^ in ^(M^) and 
respectively. We now measure the order of magnitude, as e — > 0, of each term in 
the energy: 

• Kinetic energy: the leading order term corresponds to the differentiation with 
respect to cci or X2. 



Iiv^ll 



L2 



|V/(xi,a^2)n3(£X3)pda;: 



_Q-1 



Potential energy: the leading order term corresponds to the x^ component. 



V{x)\v{x)Ydx^e'' / xi\f{xuX2)V\g{ex:i)Ydx 
Jr3 



-Q-3 



Cubic nonlinear term: 
Finally we compute: 



mi, ^2 



G{k 



?3 



G{k 



d^ 



mu^2)\'\G{^s)\'dC = o{e''^-') 



where we have used Lebesgue's Dominated Convergence Theorem. We therefore 
have: 

If we choose a < — 2, then the leading order term is the third one, and the lemma 
follows. □ 

As a consequence we are now able to prove finite time blow-up for a certain class 
of initial data. 

Theorem 5.2. Let d = 3 and G E. Denote u = minwj and assume that 

3E ui^xip\\l2. 

Then the solution ijj blows up in finite time, i.e. < oo in Corollary ^.1. More 
precisely, we can estimate ^ 7r/(2u;). 
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From Lemma 5.1 we know that for Ai < we can always choose initial data 
such that E < Q and thus enforce finite time blow-up. In other words, blow-up 
may occur even in situations where we have a defocusing (local) cubic nonlinearity 
(Ai > 0) and a positive coupling constant A2 > (the physical case), provided that 
Ai < 3-A2 holds true (and the initial energy is sufficiently small). Recalling the 
well known fact that when Ai > and A2 = 0, finite time blow-up cannot occur 
(i.e. = 00), this shows that the presence of the dipole-tcrm may indeed cause 
collapse of the Bose-Einstein condensate. 

Proof. The proof is based on the virial computation. Set 

V{t)^ f \x\^m,x)\^dx. 
Then, following ['*] (with slightly different notations), we have, since K is even: 
y{t) ^ AE + Xi\\iP{t)\\l, - 8 [ Vix)m,x)\^dx 

-2A2^^ (^(^K+^x-^K^ H^l') {t,x)m,x)\''dx. 
We again use Planchcrcl's formula with p = \t/j\'^ to write 

/ (^{k + ^x-\/K^*p^ {t,x)p{t,x)dx ^ J (^/? + ix^) \mfd^ 

and compute, for d ~ 3, 

x'^iO = - div (^K{0) = -3K{0 - ^ ■ Vif(C). 



From Lemma 2.3, we infer ^ • VA'(^) = and thus 

yit) = 4E+j^\\pm%-& [ v{x)m,x)\'dx + -^ [ Kmmi'd^- 

This can be rewritten as 
with LO = min ujj and 

fit) = AE-iY. (^1 -^') / x'^m,x)\'dx+-^ I (Ai + A2/?(0) \M)?di- 

Recalling the energy as written as in (5.1), we first note that the source term f{t) 
can be estimated via 



fit) ^^E+ -i^ / (ai + A2/?(e)) m)\'d^ 



(27r)3 

^ 6^; - II VV'(i)llL2 - 2 / V{x)\ij{t, x)\^dx ^ 6E. 



On the other hand we have 

, , , ,sin(2tjt) /■* sin(2w(t- s)) ^, , , 
2/(t) = y(0)cos(2a.t)+y(0)^-^+ ^ -f{s)ds. 

Suppose now that > 7r/(2ti>). Then we can consider t ^ 71/(201) in the above 
relation. This yields, since the sine function in the last integral remains non-negative 
on [0,^]: 

/ n\ , , r^^^-^ sin (tt - 2u}s) , , , 3E 
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We note that the left hand side must be positive, since tp{t, •) € E. By assumption, 
the right hand side is non-positive. This yields a contradiction and hence the 
result. □ 

Remark 5.3. To shed more light on the conditions for blow-up, consider the case 
of an isotropic trapping potential: ujj = to for 1 ^ j 3. The condition of 
Theorem 5.2 then reads 

^liv^ili. + + '^Mh + ^ {K, * 1^1^) {xmx)fdx ^ 0. 

In view of Plancherel's formula, this is equivalent to 

m^wh+Lo'Wxipwi. + (^^3 (ai + A2X(o) m)\'d^ ^ 0. 

From Lemma 2.3, when A2 > (the physical case), this is possible only if Ai < ^X2. 



3 

We already know that when Ai ^ ^^2; then finite time blow-up cannot occur. 
On the other hand we have just seen that blow-up occurs in the case Ai < ^A2, 
provided the initial energy is sufficiently small, say non-positive. What remains 
open therefore is the case of (large) positive initial energy in the unstable regime. 
Unfortunately we can only give a partial answer to that. 

5.2. Global existence for d = 3 (unstable regime). Here we shall show that 
global in time existence is possible in the case Ai < -^A2 under some additional 
assumptions. 

Proposition 5.4. Under the same assumptions as in Proposition 3.4, d — 3, 
suppose that in addition Ai < ^A2 and A2 ^ 0. Then there exists a Cq > 0, 
independent of Xi and A2, and a (positive) constant 

Co 



Cn = 



A2-Ai)'^ 



such that, if < E < Cq and ||V(y5j|^2 < Cq, then = 00. 

Of course, this result can only give an additional insight in situations where 
a;^||x(y9||^2 < 3Co. The picture becomes a bit clearer, though, if one ignores the 
harmonic confinement for a moment, i.e. set tOj = 0. Then we know that in the 
unstable regime Ai < 3-A2 blow-up occurs as soon as the total initial energy is 
non-positive. On the other hand, if the total initial energy is positive but not too 
large, i.e. smaller than Cq, then global in time existence still holds, provided the 
initial kinetic energy is also smaller than Cq. 

Proof. We resume the same approach as in the proof of Theorem 4.4, from which 
we now get 

wvmwh ^ 2i?+^ (j-\2 Ai) \m\'d^ - 2i?+ (^fx, Ai) wmwU- 

The Gagliardo-Nirenberg inequality then yields 

(5.2) ||v^(i)||i. ^2£; + c(^yA2-Ai) mt)\\L4^m\\h- 

On the other hand, the conservation of mass implies: 



\Vm\\h ^^E + C {^X2- Xi ] VM\\Vm\\L2- 



3 



With these two estimates in hand, the result follows from a bootstrap argument 
(see e.g. [2]): 
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Lemma 5.5 (Bootstrap argument). Let f ^ f{t) be a nonnegative continuous 
function on [0, T] such that, for every t € [0, T], 

/W ei + £2/ (^)^ 
where £1, £2 > and 6 > 1 are constants such that 



(6'£2)i/(e-i) ' ' ' " (6i£2)i/(''- 
Then, for every t G [0, T], we have 



£1. 



We apply the lemma with fit) = |j VV'(t)||i2, £1 = 2E, £2 = C (^Aa - Ai) 
and ^ = 3/2. Note that we therefore have to assume that the energy is positive. 
Otherwise, (5.2) yields no information any way, and Theorem 5.2 shows that finite 
time blow-up occurs. □ 



6. Dimension reduction 

We shall mainly follow the ideas of [h] where such an analysis has been rigorously 
performed for A2 = (no dipole effects) in the case of modulated ground state intial 
data. We also remark that the case of general initial data has been treated in the 
remarkable paper [4]. In our case, due to the presence of the non-isotropic dipole- 
kernel, we can distinguish two main cases: The reduction from d = 3 to an effective 
one-dimensional model in the dipole direction n = (0, 0, 1), and the reduction to an 
effective two-dimensional model perpendicular to n. We shall sketch the adaptation 
of the approach in [5] to the present context. It is very likely that the analysis of 
[4] can be adapted to the dipole case, but wc shall not pursue this question as it is 
beyond the scope of our work. 

6.1. Formal derivation of the one-dimensional model. Let us start with the 
first problem of deriving an effective one-dimensional model in the dipole direction. 
To this end we write the general model (3.1) in the form 



1 , , 1 / n X 1 9 3j 



(6.1) tdtij + 2 = ^ [ujlf + u^lfj ij + Lolf^ + Ai l^fip + A2 {K * i^. 

The parameter s is positive and small, its smallness modeling a strong confinement 
in the first two directions. Introduce a change of variables via 

ip{t,xi,X2,X3) = tp^ (t, yja^s) ■ 

Then equation (6.1) is equivalent to: 



1 91 1 
where H2 is the two-dimensional harmonic oscillator 



(6.2) id^^^ + i^^p^ = j-H2r + ^Iy"^' + r + ^2ic'r , 



1 f \ 2 -^1 I 2 •^2 



and /C^ is the convolution operator 

lC^{t.x)=[ (^^1 - V^f + (^^2 - v.? - 2(X3 - |^.( )|2 ^ 

U iiex^^y^r + iex2-y2r + {x3-y3rf^ 
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For £ e N, we shall denote by xe{xi,X2) the (normalized) eigenfunction of H2 corre- 
sponding to the eigenvalue G . In particular the ground state xo corresponds 
to fiQ = ^((jJi + UJ2) and is explicitly given by 



;,„(^,,^,)=^^e-(--?+--^)/^ 

IT 

As in [5], we consequently seek a solution to (6.2) in the form 

(6.3) ip''{t,xi,X2,X3) = e~^^°*^'' xoixi,X2)u{t,X3). 

This type of ansatz obviously requires well-prepared initial data (i.e. concentrated 
on the ground state of H2). A formal multiple scales expansion for il'^ as e 
then yields the consistency relation 

{H2 - iioW - 0, 

plus an evolution equation for the modulation u{t, x), given by 

(6.4) idtu + tttttM = ^l^u -f Ki\u\^u + X2 {Ki * \u\^) u, 
Here we denote by 

Ki = Ai / Xo(a:;i,X2)(ia;i(ia;2, 

the effective coupling constant for the cubic nonlincarity and by 

r a;2 _|_ _ 2a;^ 

(6.5) Ki{x3) = — i -2 -^xo{xi,X2fdxidx2. 

Jr^ {xj +xj+ xl) ' 

the corresponding effective one-dimensional dipole kernel. To apply the existence 
analysis presented in §4, it is important to check if the operator ICi : u t—f Ki * u 
is bounded on L^(IR). To do so we shall show that the one-dimensional Fourier 
transform of Ki is in i°°(M) . Recall 

Xi -f- 2xo 

K{xi,X2,X3)~ 



{xf + X2 +2^3) 



Then, denoting by J-2 the partial Fourier transform with respect to the variables 
Xi and X2 only, we obtain from Plancherel's formula on M?, that 



Ki{^3)= f e"'--'^'Ki{x3)dx3 

JR 



(2 



e ^'^^^'^ K{xi, X2, X3)xo{xi, X2)'^ dxidx2dx3 
1 f e-'-^^^T2K{i,,i2,X3)T^^{^i,i2)d^^d(,2dx3 

I") JK3 

^ / ^(Ci, 6)^6^6 / e-"^«^'^2A'(a, 6, 2^3)^x3 



i2n) 
1 

(2^ 



^2xg(a, 6)A'(ei, 6, 6) rfCid^rfCs, 



where K denotes the three-dimensional Fourier transform of the kernel K. Since xo 
is a Schwartz function, the Fourier transform of Xo is a Schwartz function, thereby in 
L^(K). Consequently, the boundedness of K on M.^ (which stems from Lemma 2.1, 
or more explicitly from Lemma 2.3) implies the boundedness of Ki on R. 

The assumptions of Proposition 3.4 are thus satisfied. From Corollary 4.2, we 
know that the effective one-dimensional model (6.4) has a unique, global solution 
in S. Moreover, the same holds true if S is replaced by L^(M), as proved in 
Theorem 4.3. 
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6.2. Formal derivation of the two-dimensional model. We proceed as before. 
Starting from 

(6.6) + = (^u;!^ + Ail^lV + X2 {K . V', 

we introduce a new change of variables by 

■>p(t, Xi,X2,X3) = ip' (t,Xi,X2, ^) . 

Equation (6.6) then becomes 

where i?i now denotes the one-dimensional harmonic oscillator, acting in the dipole 
direction 

& 



2 -^3 

3 

Denoting the corresponding ground state by Xoi^a) said proceeding analogously 
as before we arrive, instead of (6.2), at the following effective equation for the 
modulation u(t^ xi, X2)' 

(6.7) idtu + - ( + —-2 - {ujIx\ + ujIxI) u + ki\u\^u + \2K2u, 



2 \dxl 



where now 



Kl^ Xi I Xoi^3)dX3, 

is the new coupling constant and 



K2ixi,X2)= I ^ 7 ^^X0{x3ydx3 

JR (Xf + X2+ Xg) ' 

is the effective two-dimensional dipole kernel. By the same computation as above, 
we sec that T2K2 G L°°{R.'^), and thus the assumptions of Proposition 3.4 are 
satisfied to guarantee a well-posed initial value problem in E, at least locally in 
time. 

6.3. A rigorous result. For completeness we shall finally state a rigorous math- 
ematical result and sketch the corresponding proof which follows the lines of [5]. 
We only consider the case Ai ^ ^-^2 ^ (the stable regime), which allows for a 
somewhat shorter argument and again refer to [5, 4] for more general statements. 
Moreover, for notational convenience, we shall state the result only for the case of 
(6.4) (one-dimensional model in dipole direction), but exactly the same approach 
can be followed for justifying (6.7). We consider initial data for -0^ which do not 
depend on e, and refer to [•">, 4] for a discussion on this assumption. 

Theorem 6.1. Letd ~ 3, Vd be quadratic and Xi ^ ^X2 ^ 0. Denote by4>'^{t) e E 
and u{t) e E the unique solutions of (6.2) and (6.4), respectively, with 

i^^\t=o = Xo{xi,X2)uq{x3) , u\^^^ = uo{x3), Mo e E. 

Then, for any T < 00, there exists Ct > such that 



sup 

te[o,T] 



s$ Cts. 
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Proof (Sketch). The energy associated to (6.2) can be written as 



2(27r) 

where (•,•) denotes the scalar product in i^(R'^), and iJi, H2 denote the one- 
dimensional and two-dimensional harmonic oscillator operator, respectively. In the 
following we denote by 11^ the orthogonal projector onto the eigenspace correspond- 
ing to (2^3) s-iid define 

V'KijXs) = e'^"*/''' / 'ip''{t,xi,X2,X'i)xii{xi,X2)dxidx2. 



Recalling the ansatz (6.3) for il]'^{t) and using the fact that the eigenfunctions 
{xi\i&i form an orthonormal basis of L'^{W^) we note that Xo|t=o ~ '^|t=Oi since 
IIxoIIl^ = 1- Thus, we can rewrite 



{r{t),H2r{t)) = 5^(/i£ -Mo)ll#(i)lli. +A^oh(0)! 



Keeping in mind that, by assumption (i/)^(0), i?2'/''^(0)) = Mo||w(0)| 
servation of energy implies 

^ 00 

i?red = {r(t),H,rit)) + ^ Y.^^l, ~ A^o)||#(i)|li2 



the 



2(27r)3 

where the reduced initial energy is 



47r 

Ai + — A2 ( 3 



£^rcd= (^"(0),i/lV/(0)) + 



1 



2(27r) 



47r , 



^1 + — ^2 ( 3^ - 1 ) ) \m)Vd^ 



by assumption on the parameters Ai, A2. On the other hand, by exactly the same 
argument we infer from the energy conservation above that 



^ 00 



-I- 



2(27r)3 



Ai - ^A2 ) m)\'dc 



All the terms in this equation are non- negative and thus, we immediately obtain 

00 

e=i 

Denoting by Hq the projection onto the eigenspace generated by the ground state 
Xo, we can thus estimate, since V'oXo = e'''"*''^ ^oi^'^ , 



L2(R3) 



11(1 - no)7/^^(i)!l^.(R3) + e-^''o*/^\o(V'g(0 - nit)) 

= 11(1 - noWmmn^) + - • 



L2(R3) 
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For the first term on the right hand side, we have 

OG -.00 

!l(i - no)vnt)l!i.(K3) = E llV'IWIli^ < E(/^^ - Mo)IIV4(t)lli. ^ Ce'. 

Applying one derivative with respect to xi or X2, Sobolev embedding shows that 
the above estimate remains vahd if the L^(R^)-norm is replaced by L^^ (L^^ ^^), 
for 2 ^ p < 00. 

For the term -ipQ — it, more computations are necessary, and we shall emphasize 
the main arguments only. First, as in [5], we check that i/jq solves 

where 

Recall that u solves 

idtu = Hiu + Xi\u\'^u / Xoixi,X2)dxidx2 + X2 {Ki * 

where Ki is given by (6.5). Introducing the error = ipQ — u, we can write 
idtw^ = iJiw^ + + XiSf + X2S2, 

where 

SI = e'^«*/-'no (/C^V^^) - (i^i * ivgp) ^g- 

The end of the proof relies essentially on the same estimates as in the contraction 
part of the proof of Proposition 3.4 in the case d = 1, based on Strichartz estimates; 
we refer to [5] for adaptations due to the fact that the variables xi and X2 appear 
as parameters here. Roughly speaking, the term L' is treated like a linear term 
in a Gronwall lemma, even though this not so clear unless one goes through the 
details of computation (which we shall not do): for tiny time intervals, the term 
corresponding to is "absorbed by the left hand side" . Since wc consider only 
finite time intervals, we proceed this way a finite number of times, so the size of w'^ 
is dictated by the size of the source terms Sf and S*! in suitable mixed time-space 
norms. 

In view of the above remark on the estimate for (1 — no)'0'^, Sf is 0{e) in the 
spaces we need; this point is exactly the same as in [")], so we leave it out. To 
conclude, we simply consider the case of Sf. Recall 

7e3 {{exi - yi)^ + {ex2 - V2Y + (a^s - V:iY) ' 

Write 

W{t, y)|2 = |V>^(t, y)|2 _ \ii,r{t, y)P + |no^^(t, y)\\ 
Up to an error of order Oie)^ we can replace /C^ by 

r-e f (exi - j/i)^ + (£0:2 - y2)^ - 2(x3 - ya)^ 2 , 

^1= / — -. — ; -;;^\Mt^y3)xoiyi,y2)\ dy 

Jr^ {{exi - ?/i)2 + {ex2 - y2)^ + {xs - y^?) ' 

Essentially, the only new term to estimate is 
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Up to an extra "Gronwall term", we can replace this by Ilo/Cf ~ Ki* \uxo\'^, where 

^1= / ; ; -T^k(i,2/3)xo(yi, 2/2)1 dy 

Jr^ {{exi - yiY + {ex2 - 2/2)^ + (2^3 - V:iY) ' 

The difference Ilo/Cf ~ Ki* \uLp\^ =: Plt.x^) reads 
1^ = I S''{xi,X2,X3,yi,y2,y3)\uit,y3)\'^xa{yi,y2fxaixi,X2)'^dyidy2dy3dxidx2, 



where 

5" = K{yi- exi,y2 - £2:2,2/3 - 2:3) - (2/1, 2/2, 2/3 - 2:3) 

{exi - yif + {ex2 - y2f - 2(2:3 - y3f 2/? + 2/1 - 2(2:3 - y3f 



{{ex, - + {ex2 - y2? + (^3 - y3?f'^ {yl + yl + (^3 - y3Yf'^ ' 
Plancherel's formula with respect to the variables yi and y2 yields: 

^'(*'^3) = 7^ / (e--(-^''^+-^''^) - 1) ^2/^(^1,772,2/3 - X3):MxI)('7i,'72) 
X I w (t , j/3 ) p xo (2: 1 , 2:2 ) ^ dr/i fi?72 dj/3 c?2; 1 (ix2 . 



Now computing the L'^ norm of 1^ (as a function of 0:3), the boundcdncss of the 
three-dimensional Fourier transform of K and the boundcdncss of u in I? (M) yield 

[ e— (-■^''1+-=''^) - 1 V2 (Xo) {m,m)xo{xi,X2f\^ drndmdx,dx2. 

Recalling that xo G '^(M^), the standard estimate |e'^ — 1| ^ |0| for 6* € M shows 
that the norm of r is 0{e), uniformly for t ^ 0. Therefore, its {[0 , T]; L'^) 
norm (the norm that appears in energy estimates, which are a particular case of 
Strichartz estimates) is of order 0{eT), and the proposition follows. □ 
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